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ABSTRACT 

This is a short write-up of lectures at UK Theory Institute, Swansea UK, Sept. 99; and XVII 
School "QCD: Perturbative or Non-perturbative?" Lisbon, Portugal, Oct. 99. The covered topic 
0^ ' include (i) discussion of the scales at which perturbative description should be changed by the 

non-perturbative one. A number of various examples are considered, leading to different scales, 
related to instantons or confinement effects; (ii) Instanton vacuum as an Instanton Liquid; (iii) 
The phase diagram of QCD at finite temperature and density, and especially description of a 
^ , recent progress in Color Superconductivity. 
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1. Introduction 



I hardly need an explanation of why I have devoted so large portion of these lectures to such 
introductory subject as "scales" of non-perturbative QCD. First of all, it answers the question put in 
^ I the title of the Lisbon school. The second reason: somewhat surprisingly to me, naive simplistic ideas 

■ originating in 70's - the picture hadrons as some structureless perturbative "bags" , with very soft 
non-perturbative effects appearing only at its boundaries at the scale 1 f m - are still alive and well, in 
spite of multiple evidences to the contrary, pointed to existence of some substructures inside hadrons. 
The lectures in such schools is probably the right place to collect at least few simplest evidences for 

0> ■ that in one place. 

0^ I Indeed, the perturbative QCD naturally points toward the Aqcd ^ 200 MeV, the position of 

the so called Landau pole, as the momentum scale where it becomes inapplicable. Since the inverse 
Qh! of it, 1 fm, is a typical hadronic size, such cutoff indeed seemed natural. The basic approach to 

non-perturbative physics which has started in 70's, the QCD sum rules, have been based on such 
^ I a picture. If the non-perturbative fields be so "soft" , the derivative expansion, known as Operator 

• Product Expansion (OPE), should be reasonable convergent. Furthermore, smooth connection between 

hadronic and partonic approaches hoped to be possible. 

However it was realized in 80 's that the non-perturbative fields actually form structures with 
. ^ sizes significantly smaller than l-,fm. Two most important non-perturbative objects are instantons 

■ (with the average size p ^ l/S/mtl) and confining strings (with the radius of the gluoelectric field only 
about r string ~ l/5/?n, see fig. 2b ). the ratio of those to hadronic size R enters in the 4-th power for 
p and in the 2-th for r string, so in both cases only a small fraction, few percent, of hadronic volume is 
occupied by them. The simplest picture of it is that of constituent quarks (or, possibly, scalar diquarks 
as well) of the size ~ p, connected by strings. 

Can we see it in some observables? The simplest thing is to monitor the onset of the non- 
perturbative effects at small distances/ large momenta where these phenomena just become compa- 
rable to perturbative effects, and, furthermore, the OPE expansion in derivatives start to fail. Where 
it happens would define the scales we speak about, as you see the number depends on a particular 
physical problem considered: several different situation which have been identified so far will be dis- 
cussed in this section. For pedagogical reason, we do it in two rounds: at the first one (section 2) we 
present the "naive" predictions and, in many cases, see if those really work. 

In the second round (the next sections 3,4) we discuss the explanations (to the extent they are 
understood today). The conclusions we draw from these observations are that in several observables 
already the distances r — 0.1 — 0.2 used in these studies are already large enough to be outside 
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the validity domain of the original quark-gluon description, even reinforced by the OPE corrections. 
Effective approaches like Interacting Instanton Liquid Model (IILM) or Abelian Higgs Model (AHM) 
should rather be used here, and those indeed provide at least semi-quantitative explanations of the 
non-perturbative effects. 

One of the reasons the "naive" point of view mentioned above has survived for so long is that 
many applications deal with a special case of vector currents. Those are related with such long-studied 
reactions as the e~^e~ hadrons and deep inelastic scattering. The QCD sum rules work well for 
this cases, e.g. the "next twist" correction to deep inelastic scattering are surprisingly small, etc. (We 
will show some of this in section 5, when we will discuss vector correlators.) However, if one looks 
more systematically at various channels (see e.g. reviewllj), it is clear that the behavior of the vector 
channels is in fact a single exception. Cancellation of all corrections in this case is observed but not yet 
fully understood. It seems like vector currents simply cannot see these vacuum/hadronic substructure 
at all. 

We go into more detailed discussion of the instanton-induced effects in chapter 4. It is shown 
there that consistent treatment of those, to all orders in 't Hooft interaction, has been worked out. 
The so called Interacting Instanton Liquid Model (IILM) reproduces all correlation functions known 
from phenomenology (hadronic masses and coupling constants) and/or direct lattice measurements. 

We continue to discuss instanton-induced effects in the last chapter 5, but now in connection 
with "new frontiers" in QCD leading to qualitatively new phases. Among those are Quark-Gluon 
Plasma (QGP) at high T and high baryonic density; Color Superconductivity (CSC) at low T and 
high baryonic density, and the Conformal Phase at larger number of flavors. 

2. The Scales of Non-Perturbative QCD 

2.1. The Chiral Scale 

Historically the first is the so called chiral scale A^ ~ 1 GeV ^ the upper Iwiit of low energy 
effective theories such as effective chiral Lagrangians or Nambu- Jona-Lasinio modcLa Although it does 
not follow that it is also the low boundary of perturbative treatment, it cannot at least be lower. 

These effective theories appeared in 60's, before QCD was even invented. Their parameters 
were simply fitted to observations: e.g. the first Weinberg term of the chiral Lagrangian contains 2 
derivatives and the pion decay constant /^r, the next have 4 derivatives and the so called Leutwyler 
couplings, etc. The "chiral scale" appears when one asks at which momenta all terms of the chiral 
Lagrangian become of the same magnitude. 

The NJL model is somewhat more microscopic: it does not start with the Goldstone modes 
(pions), but derives instead the chiral symmetry breaking and Goldstone modes starting from some 
hypothetical 4-fermion interaction. I would not describe this model here, and only make few comments 
about it. (i) It was the first bridge between theory of superconductivity and quantum field theory: in 
a way, it first shown that the vacuum of the of strong interaction should be truly non-perturbative. 
It has a non-zero quark condensate and a gap 330-400 MeV, known as "constituent quark mass"; 
(ii) The model has basically 2 parameters: the strength of its 4-fermion interaction G and the cutoff 
A ~ .8— IGeV. The latter regulates the loops: this model is non-renormalizable. The latter parameter 
can also be seen as the inverse size of the constituent quarks; (iii) The interaction on which NJL model 
was based is hypothetical, and we should ask now if QCD actually generates it. People have been first 
trying to explain the NJL forces as being due to one gluon exchange diagrams. However it became clear 
by now that it is not so, and the true origin of the NJL-type short-range quark-antiquark attractions 
comes from the instantons; (iv) The symmetries of the resulting 4-fermion effective interaction, known 
as 't Hooft Lagrangian, are not the same as of the original NJL model. They explicitly break this 
chiral U(l) symmetry (rotating all fermions by exp^icf) * 75)), as a result of which the SU{S) / singlet 
meson rj' is not a Goldstone boson. NJL model was wrong at this point. 

*Including neutrino scattering: the leading twist operators providing moments of the structure functions are again 
vectors. The reason they correspond to the amplitude squared. Only for polarized structure functions another channel 
- the axial current - appears, and that lead to a surprise known as the "spin crisis". 
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Instead, we rather consider another incarnation of the chiral scale, defined as the energy/momentum 
scale at which the non-perturbative effects become equal to perturbative ones (in fact, the zeroth order 
in as) in = correlation functions. 

Let us consider two currents separated by space-like distance x (the spatial distance, or Eu- 
clidean time) and consider the correlation functions of the type 

K{x) ^< T{J{x)J{0)) > (1) 

with J{x) = 'ip{x)T'ip{x). The matrix T contains 75 for pseudoscalar channels, and a flavor matrices, 
if we discuss e.g. pions. For simplicity, we will restrict ourselves in what follows to 2 quark flavors. It 
means we have 4 channels: tt (P=-1, 1=1), a (P=-hl,I=0), -q (P=-1,I=0) and 5 (or ???) (P=-hl,I=l). 

In all cases at small x we expect K{x) « -^^0(2;) where the latter corresponds to just free 
propagation of light quarks. If they are massless, the correlators are Ko{x) — 12/ (tt^x^), basically the 
square of the massless quark propagator. i_. 

The first deviations due to non-perturbative effects have been calculated using the OPE,EJ in 
the context of QCD sum rules. Ignoring subtleties, one can view it as just expansion in x at small x. 
For scalar and pseudoscalar channels the result for the first correction is 

|W.l + |_<(„C)=>+,^^ (2) 

If the "gluon condensate" is made out of soft vacuum filed, they are are long range as compared to 
X and therefore all their arguments can be simply taken at the point x=0. The value of the "gluon 
condensate" appearing here was estimated previously from charmonium sum rules: 

<{gG)^ >svz--5GeV^ (3) 

Thus, the OPE suggests the scale at which the correction becomes equal to the first term: 

xoPE = (384/ < {gGf >svzfl^ ~ 1.0 fm (4) 

which seems to be completely consistent with the expectations. ^ 

However, as the Novikov, Shifman, Vainshtein and Zakharov soon noticed,H this (and other OPE 

corrections) had completely failed to describe ^1 = O"^ channels. It was very puzzling, because 

for vector and axial channels it did a good job.liil That is why they asked an important question (the 

title of their paper) : "Are all hadrons alike?" . 

How do we know that this scale is wrong? Let us just do the simplest thing: consider the pion 

contribution to the correlator in question 

The coupling constant is defined as A^r =< 0|J(0)|7r > and the rest is nothing else as the scalar 
massless propagator. (We can ignore the pion mass at distances in question.) 

Because both the pion term and the gluon condensate correction are 1/x'^, let us compare the 
coefficients. Ideal matching would mean they are about the same^: 

< jgG)^ >svz 
87r2 



^TT ^ (,0j 



The r.h.s. is about 0.0063 GeV^. However (unlike the better known coupling to the axial current f.^), 
the is surprisingly large^. The l.h.s. of this relation is actually Aj = (ASGeV)'^ = 0.053 GeT^^, 

tin fact there are states other the pion, and all of the contribute positively to the correlator. So one should rather expect 
the pion contribution be somewhat smaller than the correction term, to be corrected by the contribution of other states. 
tThe reason for that is the the pion is rather compact and also the shape of the wave function is concentrates at its 
center, so that its value at r=0 is large. We return to this point in the discussion of the "instanton liquid" model. 
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or about 10 times larger than the r.h.s. It means some other and much larger effect should explain 
deviation from the perturbative behavior. Why it has not showed up in the vector channels? We return 
to these questions in the next section. 

2.2. The Scale of the — glueball correlation functions 

It is defined in the same way as the first one, only for scalar and pseudoscalar glueball correlators. 
Instead of the quark currents we have now operators G^^ and G^j^G^^ (tilde means dual field in 
electric- magnetic sense, or convolution with £^^^5/2. The analogous OPE expressions for these two 
channelsu now starts with the next gluonic "condensate" 



Kj^M^) = ±^[1 ± ^ < g'f^'^G^.AG^,. > +...] (7) 



Again, the first term is the free propagation of gluons, it is the largest one at small x. Note that the 
power of X in both terms simply follows from dimension. Demanding that these two terms are equal, 
we get the corresponding "OPE scale" . ^ 

However, as emphasized in this paper J3 other theoretical arguments suggested that in fact the 
non-perturbative physics should start at much smaller scale. One of them is based on low energy 
thcorcn: 



/ 



d'xK+{x) = < (gG)' > (8) 



where b= ^ " ^ « 9 is the first coefficient of the beta function 



One way to see what it means is to assume that the major deviation at small distances is due 
to the contribution of the scalar glueball^ Then the integral can be related to the coupling constant 
Ao+ =< 0|G^|0+ glueball > and the glueball mass mo+ in a simple case 

^ = ^ < (gG? > (9) 



'"-0+ 

It is known from lattice and phenomenology that too+ ~ 1.6GeT^, so we now can obtain the coupling. 
Then we can go back to K^{x) at small x and see where the glueball contribution becomes equal to 
the perturbative part, to define the scale. 

Like in the pion case, both the OPE correction and the coupling produce 1/x'^ effects at small 
distances, so we again compare the coefficients. The OPE-hadron matching would work if they are 
close 

Ag+ « 8 < gV^'^G^l^G^G^, > (10) 

but in fact the value we infer from the low energy theorem gives 1.8 GeV^ for the l.h.s., which 
again significantly exceeds the l.h.s. (when the phenomenological value of the VEV of G^ operator is 
substituted.). 

More strict wayEI to say the same thing is to go to momentum representation (Fourier transform 
(g) of the correlation function. The total integral over x corresponds to the zero momentum value. 
Since we know it, we can subtract it in the dispersion relation 



SThe analogous integral for the pseudoscalar correlator is the so called topological susceptibility. 

Hxh ere are evidences that the scalar glueball is even more compact and stronger bound object than the pion, so it is a 
reasonable assumption. 

I Lattice and instanton liquid calculation of the coupling leads to a number consistent with this estimate and the low 
energy theorem is of course exactly satisfied. 
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where — —q^ > 0. The subtraction improves the convergence of the integral at large Q of the 
dispersion relation, namely 



K+{Q)_K+{0) , 1 f ^^ Imk+{s) ^^2) 



Q2 s{s + Q 

This argument shows that in momentum space the effective subtraction term becomes equal to 
perturbative one at the momentum scale ^Q±giuebaUs ^ (3 — 4) GeV, the largest of the effective scales 
of the non-perturbative QCDp*] Other set of arguments presented in the same paper have shown, that 
the same scale should appear in the pseudoscalar glueball channel as well. Why this happens in this 
case at much larger scale than for any quark channels? 

2.3. Heavy Quarkonia and Small Size Instantons as Dynamical Dipoles 

Another way of looking at the onset of non-perturbative physics is to separate color charges 
by some small distance, and see what happens. (The correlation functions considered above deal with 
two colorless operators, also separated by a small distance x.) We will call such cases generically "a 
small-size color dipoles" . There would be three different kinds of those. 

Historically the first example are states of heavy quarkonia. The non-perturbative correction to 
their energies was calculated by Voloshin and Leutwyleiiiy by OPE: SE ^< 0\g^G^i,\0 > t where 
the spatial size r ~ l/{asM) and the rotational time t ~ l/(a^M), both small for large quark mass 
M. (I am not going to bother you with exact expression, in which r appears as some dipole matrix 
element and r as energy denominators, as well by questions how small is as and what is its argument.) 

The question is whether this correction is indeed what is the case in reality. The answer probably 
exists or can be worked out using the spectroscopic information about Upsilons at hand: unfortunately 
I am not aware of any precision studies of it. 

Let me therefore move on to another kind of "dynamical dipoles", the instantons, now 4- 
dimensionally symmetric, with r ^ t ^ p. One of the reason is I know it better, but there is also 
a more scientific explanation of why one should do it instead. Instantons in general are much more 
sensitive tool, because the tunneling probability contains the perturbative/semi-classic term (as well 
as all non-perturbative corrections) in the exprment. 

(A side remark: as noticed in my paper ,113 it is in particular should be very important for fixing 
the value of A^qcd- Experimentalists made tremendous efforts to get it from scaling violation and 
similar effects containing \og{AQC'D), but still get poor accuracy. Lattice practitioners use hadronic 
masses and especially splittings of some quarkonium levels for this purpose, which are effects const * 
Aqcd: a-iid improved the accuracy. Butj-thc the density of small size instantons, calculated semi- 
classically in a seminal paper by 't Hooftjlil contains 0{Aqcd)'' where b = (117Vc/3 - 2iV//3) - 9. So 
the instanton-based determination of Aqcd should be potentially 10 jtimes more accurate!) 

Similar to Voloshin-Leutwyler correction, the OPE-based resulttil predicts the following correc- 
tion to the density of instantons of size p: 

4 4 

dn{p) = dnp,rt{p){l + ^ < Ols'G^jO > +...) (13) 

It is also second order dipole approximation: the instanton dipole is 0{p^/g). Note the generic 4-th 
power of p is also determined by the dimension of the lowest gauge invariant local operator. Note 
also the sign: it is nothing else but a generic attraction resulting from any second order perturbation. 
However both these conclusions happen to be in apparent conflict with the lattice data (see below). 

eFig.[l](a) we show recent lattice data for the instanton size distribution in pure SU(3) gauge 
(There are others but, this work includes such refinements as improved lattice action and 

** Let me emphasize it once again: it is not the glueball mass scale, but rather the so called gluon-glueball duality 
energy. Above it the spectral density can be approximated by the partonic one, below it cannot. 



5 



back extrapolation to zero smoothening.) One can clearly see, that a rapid rise at small p turns into 
a strong suppression. The former behavior is consistent with the semi-classical one-loop resultEj: 




Fig. 1. (a) The instanton density dn/dpd*'z, [fm"^] versus its size p [fm]. (b) The combination p~^dn/dpd'^^ 
in which the main one-loop behavior drops out for Nc ~ 3, Nf = 0. The points are hom the lattice workJ13 
for this theory, with /3=5.85 (diamonds), 6.0 (squares) and 6.1 (circles). Their comparison should demonstrate 
that results are rather lattice-independent. The line corresponds to the proposed expression ~ exp{—2-Kap ), 
see text. 

Sharp maximum seen in Fig.[l](a) appears at rather small < p >« 1/3 /m, much smaller that 
their spacing Rk, \ fm. This results in a non-trivial "vacuum diluteness" parameterEJ to which we turn 
in section 5, but now we are not interested in a typical instanton size but rather in their suppression. 
Therefore we re-plot the same data in Fig. [1] (b) , with the leading semi-classical behavior taken away. 
Now the maximum is gone and one finds the same suppression pattern at both sides of the maximum. 
The OPE prediction (|3|) is not seen: probably it is only true at very small p. The suppression effect 
is clearly 0{p^), and not just for small p but in the whole region. 

2.4- Small Static Dipoles 

Another way to look at these issues is to use small static color dipoles, or short strings. As time 



scale T is now unlimited, there is no OPE prediction like (13). However using just the second order 
dipole approximationcH one gets a similar prediction of the non-perturbative corrections to the static 
potential 



y(r) 



4Q;s(r) 1 



+ r^ J dT e^-'-^^^ < 0|G^.(t)C/,G^,(0)C/+|0 >) (15) 

where the field strengths are separated by the time delay t, with Ut being the appropriate parallel 
transports. Note that it is 0{r^). 
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Fig. 2. (a) The potential at small distances from,Ej the quantity ro ~ 0.5 fm denotes the so called Sommer 
scale is used in lattice works, (b) The fits of the transverse distribution of aiectric and magnetic fields in the 
string, obtained in lattice calculation, to ANO vortex solution of the AHM.u It leads to Higgs mass 1.3 GeV. 



However recent lattice data on V(r) at small have found a clear 0(r) effect instead (suggested 
previously intHi): 

Vir)=-'-^'-+aor + ... (16) 
6 r 

The small-distance tension is larger than the asymptotic one ctq ~ (4 — 5)croo- (I however had to 
warn the reader that it has rather uncertain error, because subtraction of the perturbative potential 
depends on technical details.) 

This 0(r) potential for short strings makes sense to me mostly becauS|e-|tlLere are many other 
evidences that the QCD confining strings are surprisingly thin. Lattice studiestJ'El (the latter is shown 
in Fig. 2b) show in particular that the so called "energy radius" (at which it decreases by 1/e) is about 
^i/e « .18 fm, while that for the action distribution is about twice larger. 

This observation also has many phenomenological consequences. One is just another argument 
explaining weak string-string interactions known from Regge phenomenology. Another is "hadron 
diluteness" : color field inside hadrons occupy only few percent of the volume 

{S.^jRHr - (1/5)' (17) 

This is contrary to the MIT bag model which views the whole hadronic interior to be in the per- 
turbative phase. In other words, the value for the bag constant Bmit ^ 50 MeV/ fm^ was hugely 
under-estimatedi it is i? ~ 1000 MeV/ frn? or more. (Similar but different argument was made two 
decades ago in. a) 

3. Instanton-induced effects at small distances 

3.1. The Renormalized Charge 

We have not listed this issue in the list of observations presented above, because the effective 
charge itself is not so easily observable. Nevertheless, all calculations deal with it in some way or 
another, and so it is impossible to avoid it in a theory discussion. So let us ask what constitutes 
the first non-perturbative correction to the effective charge. It is defined here as the coefficient in 
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an effective action for some external weak background field GJ^^^, in which both perturbative loop 
corrections and the non-perturbative effects are to be included. 

A very simple non-perturbative effect was suggested long ago by Callan, Dashen and Gross 
(CDG)l3: instantons are simply dipoles and therefore they are polarized by external field and produce 
some dielectric constant, exactly like atoms do. The external field is supposed to be normalized at 
some normalization scale /z, and CDG has proposed to include all instantons with size p < pmax = 1//^- 
The effective charge is then defined as: 

- = 6 In -t— ^ — 5 -/ dn{p)p\ , , . f (18) 

where b = llNc/3 — 2Nf/3 is the usual one-loop coefficient of the beta function, and dn{p) is the 
distribution of instantons (and anti-instantons) over size. 

CDG could not calculate it then, as the instanton size distribution was unknown. Plugging in 
lattice data (such as shown in Fig. la) for SU(2) gluodynamics, SU(3) gluodynamics we getEl curves 
shown in Fig.(|^). The curve for QCD with fermions is not from lattice: we do not have it yet, so 
it corresponds to IILM. Note a rapid departure from the perturbative logarithmic running coupling 
constant. 




Fig. 3. The effective charge b g'^j f{p) /Sn^ (b is the coefficient of the one-loop beta function) versus normal- 
ization scale p (ip-.units of its value at which the one-loop charge blows up). The thick solid line correspond 
to exact solutiontj for the N=2 SYM, the thick dashed line shows the one-instanton correction. Lines with 
symbols (as indicated on figure) stand for N=0 QCD-like theories, SU(2) and SU(3) pure gauge ones and QCD 
itself. Thin long-dashed and short-dashed lines are one and two- loop results. 

Is this indeed what is happening? In order to answer such questions above I have used either 
lattice data or hadroniajphenomenology: let me now use another tool from the theory tool-box, namely 
a solvable model. So wcEI are going to compare the charge in QCD to that for the N=2 superaj^mmetric 
theory, for which the famous exact effective Lagrangian was found by Seiberg and Witten.l 



E4 



877 _ A-(7T~fc^) 
9^{u) K{k) ^ ' 
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where K is elliptic integral and the argument 



e . '-v;'"'-''^'" (20) 

(w + _4A4)) ^ ^ 

is a function of gauge invariant vacuum expectation of squared scalar field 

. = -<0^>^- + ^ + ... (21) 
and a is just its VEV. For large a there is a weak coupling expansion which includes instanton effects 



Svr 2 / /2a2\ SA^ , 

log TT -^ + --- (22) 



The exact coupling blows up at u = 2A^, which means that the factor between the exact strong 
interaction scale and the perturbative one is in this theory A^o = 2'^/^A. Actually this is the ratio of 
the scale ^/u to the scale a at which the perturbatively evolved coupling (one-loop) blows up. If one 
were to account for the next term in the expansion of u, the ratio of scales is reduced to \/2\/2 + \/2. 
The fact that instanton effects can be important at such a high scale was anticipated in Ref.B and is 
presumably due to the significance of the pre-factor in instanton calculations.) 

The behavior is shown in Fig.(^, where we have included both a curve which shows the full 
coupling (thick solid line) , as well as a curve which illustrates only the one-instanton correction (thick 
dashed one). The units on both axis are chosen so that the naive one- loop running looks the same in 
both theories. And - surprise, surprise - the instanton-induced corrections seem to be the nearly the 
same as well. 

I think this comparison teaches us few lessons. The main of them: if the instanton-induced 
corrections to the charge itself becomes large, the perturbative expansion in g have to be abandoned. 
I do not know what exactly one can make out of (nearly ideal) numerical matching of these two 
theories, QCD and N=2 SUSY gluodynamics. Note however the very rapid change of the coupling 
induced by instantons, although different in formulae, but so similar in the curves. It is also of interest 
that the full multi-instanton sum makes the rise in the coupling even more radical than with only 
the one-instanton correction incorporated. It is also interesting to observe that at the scale where the 
true coupling blows up, the perturbatively evolved coupling is still not very large. Individually, the 
perturbative log and instanton corrections are well defined at this region: however they cancel each 
other in the inverse charge. This is encouraging from the point of view of developing a consistent 
expansion for the instanton corrections. The rapid rise in the coupling is also encouraging in that 
it ensures that perturbation theory is valid almost to the point where it blows up. For a consistent 
picture of QCD, in which perturbation theory still appears to be applicable at the c-quark scale, while 
the theory is non-perturbative at 1 GeV, such a dramatic effect is essential. 

3.2. The Chiral Scale in Spin-0 Light-Quark correlators 

The physical origin of the chiral scale has been traced down to instantons-induced effects^ for 
recent detailed review see.Q The main points are: (i) the instanton-induced interaction strong and 
short-range, (ii) it is attractive for the pion channel, repulsive for the 77', and (iii) (to first order) it 
does not affect the vector channels such as p, lo. It explains hath successes and failures of the OPE sum 
rules, reproducing the phenomenology also for all spin-0 correlation functionsJl3 So far, we have only 
mentioned the pion contribution, and noted that it deviates from perturbative behavior at unusually 
small distances. 

After this prelude we return to evaluation of the short distance behavior of correlation functions 
in the single- instanton approximation (SIA) .E3 The main idea is that if a; — y is small compared to the 

ttit should be noted that the first terms in this expansion have been exphcitly verified in instanton calculations.^ 



9 



typical instanton separation R, we expect that the contribution from the instanton I — closest to 
the points x and y will dominate over all others. For quark propagator in the subspace of the instanton 
zero modes this implies 

s{x,v) = Y^M^) (Tf^) i'liy) ^ (TfT^) 4Ay) ^ ^''^''^f'^y^ (23) 

\T + imJ j j \T + im J j^j^ m* 

where we have approximated the diagonal matrix element by its average, {T + im)j\ ~ N^^ + 
im)Jj, and introduced the effective mass m* {m*)~^ = ^ A~^. In the following we will use the 
mean field estimate m* = 7r/9(2n/3)^/^. As a result, the propagator in the SIA looks like the zero mode 
propagator of a single instanton, but for a particle with an effective mass m* . 

The TT and rj' correlators receive zero modes contributions. In the single instanton approximation, 
we findta 

where — x^ j^x^ + 4p^). There is also a non-zero mode contribution to these correlation functions 
which is not very important. 

Putting in standard instanton parameters, we find that at very short distances a; ~ .3 fm, this 
contributions to the correlator coming from the instanton zero modes are larger than the contributions 
included in the OPE (see Fig. 4a). 

For the 77', the instanton contribution is strongly repulsive. The flip of the sign is easy to 
understand: the instanton corrections comes only from flavor-changing diagram uu — > dd, which has 
different sign in correlators of J^ro — (U75U — (?75C?) and J^,/ — {urj^u -I- ^750?) This is because t 
instanton accounts for the \J(\)a anomaly. The SIA was extended to the full pseudoscalar nonet in. 
It was shown that simply replacing m* — > m* -t-mg gives a good description of SIJ(k>) flavor symmetry 
breaking and the tt, if, 77 correlation functions. 

In the vector p, lo correlators the zero modes cannot contribute since the chiralities do not 
match. Non-vanishing contributions come from the non-zero mode propagator and from interference 
between the zero mode part and the leading mass correction 

nf^(x,y) = Tr[7^5"-(a;,y)7M5"^(2/,x)]+2Tr[7^V'o(a;)7^J(2/)7MA(y,x)] (25) 

The latter term survives even in the chiral limit, because the factor to in the mass correction is 
cancelled by the 1/to from the zero mode. Also note that the result corresponds to the standard 
DIGA,__sp multi-instanton effects are not included. After averaging over the instanton coordinates, we 
find " 



nr(x) = nO-,/dp.(p)i|^^|-^iogi^^ (26) 



Similar to the OPE, this correction is only weakly attractive at intermediate distances. 

An interesting observation is the fact that it is the only singular term in the instanton contri- 
bution. In_feict, the OPE of any mesonic correlator in any self dual field contains only dimension 4 
operatorsEd This means that for all higher order operators either the Wilson coefficient vanishes (as 
it does, for example, for the triple gluon condensate (ff'^/°'"^G^,yG(^pGp^)) or the matrix elements of 
various operators of the same dimension cancel each other. This is a very remarkable result, because 
it helps to explain the success of QCD sum rules based on the OPE in a number of channels. In the 
instanton model, the gluon fields are very inhomogeneous, so one would expect that the OPE fails 
for X <C jO. The Dubovikov-Smilga result shows how many observables can be absolutely blind to very 
strong gauge fields, as long as they are locally self dual. 

M There is a mistake by an overall factor 3/2 in the original work. 
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3.3. Instantons and the duality scale for J — O glueballs 

Let me first explain qualitatively why in this case the scale should be larger than in the light 
quark channels. The reason is the gluonic field strength inside instantons is very large. 

The magnitude of the color fields inside the instanton is 0{l/g), and in the correlation functions 
in question (such as < G"^ {x)G'^ {y) >) it enters in the 4-th power. Therefore, the instanton correction 
is proportional to the action of the instanton squared, (Stt^/^^)^ ~ 100 where for estimate we have 
used the action to be of the order of 10. The corresponding correction with light quarks have zero 
modes instead of gluonic field strength, and its integral over space is not the action but just one, 
from the normalization of fermionic mode. So we speak about enhancement of roughly 2 orders of 
magnitude. 

For one-instanton calculation we can start from expansion of the gluon operators around the 
classical fields. In thc_Lowest order one should simply insert the classical field of an instanton into the 
operators, it leads tola 

T^SiA, ^ f, ,,8192^2 Q3 (ffiO-Q^^ 3, 1 + eM 



where ^ is defined as in (24). There is no classical contribution in the tensor channel, since the stress 
tensor in the self-dual field of an instanton is zero. Note that the perturbative contribution in the 
scalar and pseudoscalar channels have opposite sign, while the classical contribution has the same 
sign. To first order in the instanton density, we therefore find all three scenarios : attraction in the 
scalar channel, repulsion in the pseudoscalar and no effect in the tensor channel. 

The single-instanton prediction can be compared with the OPE, and it is much larger indeed. 
Furthermore, if one uses the standard p, n of the instanton vacuum, it matches it the scale inferred 
from low energy theorem very well. 

!_. Numerical calculations of glueball correlators in different instanton ensembles were performed 
in.E3 At short distances, the results are consistent with the single instanton approximation. At larger 
distances, the scalar correlator is modified due to the presence of the gluon condensate. This means 
that (like the a meson), the correlator has to be subtracted and the determination of the mass is 
difficult. In the pure gauge theory we find mQ++ = 1.5 GeV. While the mass is consistent with QCD 
sum rule predictions, the coupling is much larger than expected from calculations that do not enforce 
the low energy theorem. 

Incj we also measured glueball wave functions. The most important result is that the scalar 
glueball is indeed small, ro++ = .2 fm and determined by the size of an instanton, while the tensor is 
much bigger, sensitive to interactions at the size determined by the average distance between instantons 
1 fm. 

In the pseudoscalar channel the correlator is very repulsive and there is no clear indication of 
a glueball state. In the full theory (with quarks) the correlator is modified due to topological charge 
screening. The non-perturbative correction changes sign and a light (on the scale of the glueballs) state, 
the r]' appears. Non-perturbative corrections in the tensor channel are tiny. Isolated instantons and 
anti-instantons have a vanishing energy momentum tensor, so the result is entirely due to interactions. 

4. Confinement effects at small distances 

4.1. Instanton Suppression 

The central idea of this Section is that 0{p^) suppression of instantons can be due to a "dual 
superconductivity" ,l3 a scenario in which some composite objects condense, forming the non-zero 
vacuum expectation value (VEV) of the magnetically charged scalar field cj). My first (naive) argument 
was that in such theory, unlike the QCD itself, at leaai_there is the dimension-2 operator 

The composites may be magnetic monopoles,E^O or P-vortices, or something else: anyway one 
is lead to an incarnation of the old Landau- Ginzburg effective theory, Abelian Higgs Model (AHM), 
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describing interaction of a "dual photon" and "dual Higgs" fields. AHM was applied to the description 
of the QCD strings, as Abrikosov-Nielsson-Olesen vortices.EEl 

Before we go into details, let us point out a striking similarity between these two problems. A 
vortex is the 2d topologically non-trivial configuration, in which (j) vanishes at the center, the Dirac 
string where the dual potential must be singular. An instanton problem is in a way the previous one 
squared. The 4d picture of the fields is like two string cross sections in two orthogonal 2d planes. 
Higgs field again vanishes at the center, because in the singular gauge (the only one good for multi- 
instanton configurations) the gauge field is A^{x)'^ ~ l/x^ at the origin, acting as a centrifugal barrier. 
Since "melting" of the dual superconductor at the center is not a small modification, one generally 
cannot expect the OPE-type calculations to hold. In both problems one has first to solve for the 
field and then calculate tlii3 energy or action. Fortunately, for instantons in a Higgsed vacuum it was 
already done by 't HooftllJ: for fundamentally charged Higgs the answer is 

AS* = 47r2p^| < ^ > 1^ (28) 

Note that it leads to the 0{p^) suppression law we need to explain Fig.[l](b), and that AS should not 
necessarily be small. We return to speculations on the exact nature of the Higgs and the dual photon 
fields of the Landau- Ginzburg model (needed to evaluate the strength of the effect) below. 

Encouraged by this, we return to instantons and try to apply the same reasoning. Since both 
the 't Hooft correction (^) and the string tension (^ij) scales as the Higgs VEV squared, we expect 
qualitatively that 

— = —\pertexp{-C(Tp ) (29) 

where C is some numerical constant. In order to find C one has to identify the scalar and the dual 
photon fields of the AHM, and explain how-they are coupled to the colored gauge field of the instanton. 
In the AHM treatment of the QCD string^ the magnetic field of the dual photon is identified directly 
with the color-electric gauge field inside the string. It does not create problems because this electric 
field can be considered Abelian. 

Applying the same ideas for instantons, let us first note that their self — duality helps: because 
electric and magnetic fields are identical, the "magnetic" potential and the original one are 
also the same. However both are intrinsically non-Abelian, so only a particular component (or a 
combination of those) can be identified with the Abelian "dual photon" of the effective theory. In 
other words, an Abelian projection is inevitable, and there is no unique or preferred way to do it. 
Lacking better ideas, we simply do what lattice people do: just select one of possible projections 
and see what happens. Clearly, selecting Higgs field interacting with a particular component of the 
gauge potential means breaking the gauge group (which the vacuum of the Standard model does and 
that of QCD does not). But we proceed anyway, simply re-scaling the dual fields in a way that their 
Lagrangian (1) matches the 't Hooft one. If we do so, it leads to identification < $ {2/3)<j>Q, or 
the constant C in (|2^) to be C = 27r. Putting it all together, we can now compare^ this result to the 
exponential suppression. The corresponding curves are shown in Fig.l, and it works very well. 

Brief summary. The "dual superconductivity" leading to confinement in the QCD vacuum 
seems to be surprisingly robust. Instanton suppression and/or string tension considered before both 
fix the AHS Higgs VEV Tathei accurately. The string size and small-r potential provide hints that the 
Higgs and dual photon masses are large, in the 1-1.5 GeV range. Their exact nature remains unclear^. 

Outlook: one may test our suggestions by comparing instanton suppression in theories with 
variable number of colors and/or flavors. Unfortunately available data for the SU(2) color group or 
SU(3) with dynamical fermions are not yet good enough to do so. 

* The standard string tension value <j = (AbOMeV)'^ is traditionally used to set absolute scale of the lattice data for 
non-physical theories like for the pure SU(3) gauge theory simulations we used: so the instanton sizes expressed in fm 
is consistent with it. 

* So the status of AHM Higgs is, ironically, not that different from Higgs particle of the Standard Model. 
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Another challenging set of questions is related with instanton suppression at non-zero tem- 
peratures T and/or densities. At high T or density, in the quark-gluon plasma phase, the answer is 
clear: the instanton electric fields are again suppressed, but now by the usual Debye screeninga It 
leads to a factor e*p(— a(T, ii)p^) similar to the one discussed above, where the coefficient a(T, /i = 0) 
was calculated inEa and then generalized to the fJ- ^ in.lll Note that confinement is not completely 
gone: it remains for spatial Wilson loops. The most interesting point is what happens close to the 
deconfinement transition. Since 2 and 3 color gouge theories have it of the second and the first kind, 
respectively, a detailed study of the instanton suppression at T Tc is of great interest. 

4-2. Long and Short Strings 

Now we briefiy review applications o£-the dual superconductivity idea to the QCD string, or the 
ANO vortex line, done in a series of papers .Ea Among clear successes of this approach is: (i) prediction 
of weak string-string interaction, putting it around the boundary of type I and II superconductivity; 
(ii) prediction of a whale set of potentials other than central. Both agree well with available lattice 
data, for a review see.t3 |_. 

The effective Lagrangian used inE3 is 



L ^ 7;hid^C^ - 9.C^f+ (30) 



-[- 



\2\2] 



where we have omitted interaction with quarks at the ends. C\ is dual color potential coupled to Higgs 
with magnetic coupling gm — 27r/(7. Assuming that we are exactly at the boundary of the type I and 
II superconductivity, the masses of the Higgs and the "dual photon" are equal — Mc — gm4>o- 
The (classical) string tension is directly related to Higgs VEV 

<j = Y^M' (31) 

In effective dual modelil the string width is related to masses of dual photon and Higgs, being 
the large non-perturbative scale of the 3-ed kind we are speaking about. The data mentioned put it 
in the "glueball mass range", around 1.3 GeV according to Fig. 2b (It is difficult for me to access the 
error involved.) 

5. QCD Vacuum as the Instanton Liquid 

5.1. The qualitative picture 

It was already mentioned that, as pointed out in,0 the typical instanton size is significantly 
smaller than their separation R = n^^l^ l/^i, or 

1/3 fm. (32) 

If so, the following qualitative picture of the QCD vacuum emerges 

1. Since the instanton size is significantly smaller than the typical separation R between instantons, 
p/R ^ 1/3, the vacuum is fairly dilute. The fraction of spacetime occupied by strong fields is 
only a few per cent. 

2. The fields inside the instanton are very strong G^^ 3> ^"qcd- This means that the semi-classical 
approximation is valid, and the typical action is large 

5o = Sn^/g'^ip) - 10 - 15 > 1. (33) 

Higher order corrections are proportional to 1/5*0 f^nd presumably small. 
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3. Instantons retain their individuality and are not destroyed by interactions. From the dipole 
formula, one can estimate 

\5Sint\^{2-3)<^So. (34) 

4. Nevertheless, interactions are important for the structure of the instanton ensemble, since 

expjJS'intl ~ 20 > 1. (35) 

This implies that interactions have a significant effect on correlations among instantons, the 
instanton ensemble in QCD is not a dilute gas, but an interacting liquid. 

The aspects of the QCD vacuum for which the instantons arc most important arc those related 
with light fermions. Their importance in the context of chiral symmetry breaking is related to the 
fact that the Dirac operator has a chiral zero mode in the field of an instantons. These zero modes 
are localized quark states around instantons, like atomic states of electrons around nuclei. At finite 
density of the instantons those states can become collective, like atomic states in metals. The resulting 
de-localized state corresponds to the wave function of the quark condensate. 

5.2. How to calculate? 

In this lectures there is no place for technical details, which may be found in the original papers. 

Rather wc try to explain the logical structure and possible alternatives of the theory. For that reason, 
discussion of possible methods to calculate something with instantons is put into the introduction, 
with only results of the calculation reported below. 

In the QCD partition function there are two types of fields, gluons and quarks, and so the first 
question one addresses is which integral to take first. 

(i) One way is to eliminate gluonic degrees of freedom first. Physical motivation for that may 
be that gluonic states are heavy and an effective fermionic theory should be better suited to derive an 
effective low-energy theory. In the instanton framework, this approach was started by 't Hooft who 
discovered that instantons lead to new effective interactions between light quarks. We will present the 
explicit form of such four-fcrmion interaction for two-flavor QCD below. 

It is a well-trotted path and one can follow it along the development for a similar four-fermion 
theory, the NJL model. One can do simple mean field or random field approximation (RPA) diagrams, 
and find the mean condensate and properties of the mesons. Unfortunately, it is difficult to do more. 
For example, baryons are states with three quarks, and using quasi-local four-fermion Lagrangians 
for the three body problem is technically a very difficult (although solvable) quantum mechanical 
problem. There were no attempts to sum more complicated diagrams. 

(ii) The opposite strategy can be to do fermion integral first. It is a simple step, because they 
only enter quadratically, leading to a fermionic determinant. This is the way lattice people proceed. 
In the instanton approximation, it leads to the Interacting Instanton Liquid Model, defined by the 
following partition function: 

^= ^ WInJ. n [dnid{pi)]eM-Sint)l[det{D + mf), (36) 

describing a system of pseudo-particles interacting via the bosonic action and the fermionic determi- 
nant. Here dfli = dUi d'^Zi dpi is the measure in color orientation, position and size, associated with 
single instantons and d{p) is the single instanton density d{p) = drij j/dpdz. 

The gauge interaction between instantons is approximated by a sum of pure two-body in- 
teraction Sint = ^'^i^j Sinti^ij). Genuine three body effects in the instanton interaction are not 
important as long as the ensemble is reasonably dilute. Implementation of this part of the interaction 
(quenched simulation) is quite analogous to usual statistical ensembles made of atoms. 
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As already mentioned, quark exchanges between instantons are included in the fermionic deter- 
minant. Finding a diagonal set of fermionic eigenstat es of the Dirac operator is similar to what people 
are doing, e.g., in quantum chemistry when electron states for molecules are calculated. The difficulty 
of our problem is however much higher, because this set of fermionic states should be determined for 
all configurations which appear during the Monte-Carlo process. 

If the set of fermionic states is however limited to the subspace of instanton zero modes, the 
problem becomes tractable numerically. Typical calculations in the IILM involved up to N~ 100 
instantons (+anti- instantons): it means that the determinants of N x N matrices are involved. Such 
determinants can be evaluated by the ordinary workstation (and even PC these days) so quickly, that 
straightforward Monte Carlo simulation of IILM is possible in a matter of minutes. On the other hand, 
expanding the determinant in a sum of products of matrix elements, one can easily identify the sum 
of all closed loop diagrams up to order N in the 't Hooft interaction. Thus, in this way we take care 
of (practically) all orders! 

5.3. Correlation functions 

There is no place here to explain the details of how we performed calculations in the Interacting 
Instanton Liquid Model (IILM): they are done by Monte Carlo, in a way not so different from standard 
methods used in other statmech models. Different versions of the model (mentioned in figures below 
as IILM(rat) etc) differ by a particular ansatz for gauge field used, from which the interaction is 
calculated. Note also, that these figures mostly contain also a curve marked "phen" : this is how the 
correlator actually look like, according to phenomenology. 




0.5 1 1.5 

T[fm] 



Fig. 4. Pion correlation function in various approximations and instanton ensembles. In the top figure we show 
the phenomenological expectation (solid), the OPE (dashed), the single instanton (dash-dotted) and mean field 
approximations (dashed) as well as data in the random instanton ensemble. In the bottom figure we compare 
different instanton ensembles, random (open squares), quenched (circles) and interacting (streamline: soUd 
squares, ratio ansatz solid triangles). 
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Fig. 5. Eta prime meson correlation functions. The various curves and data sets are labeled as in Fig. 



We simply jump to show a sample of results. Correlation functions in the different instanton 
ensembles were calculated in several papers (see original refs inQ). Some of them (like vector and 
axial-vector ones) turned out to be easy: nearly any variant of the instanton model can reproduce well 
the (experimentally known!) correlators. Some of them are sensitive to details of the model very much: 
two such cases are shown in Figs. ^|]. The pion correlation functions in the different ensembles are 
qualitatively very similar. The differences are mostly due to different values of the quark condensate 
(and the physical quark mass) in the different ensembles. Using the Gell-Mann, Oaks, Renner relation, 
one can extrapolate the pion mass to the physical value of the quark masses. The results are consistent 
with the experimental value in the streamline ensemble (both quenched and unquenched), but clearly 
too small in the ratio ansatz ensemble. This is a reflection of the fact that the ratio ansatz ensemble 
is not sufficiently dilute. 

In Fig. H we also show the results in the p channel. The p meson correlator is not affected by 
instanton zero modes to first order in the instanton density. The results in the different ensembles are 
fairly similar to each other and all fall somewhat short of the phenomenological result at intermediate 
distances x ~ 1 fm. We have determine the p meson mass and coupling constant from a fit. The 
p meson mass is somewhat too heavy in the random and quenched ensembles, but in very good 
agreement with the experimental value mp = 770 MeV in the interacting ensemble. 

Since there are no interactions in the p meson channel in the RPA, it is important to study 
whether the instanton model provides any binding at all. In the instanton model, there is no confine- 
ment, and nip is close to the two (constituent) quark threshold. In QCD, the p meson is also not a 
true bound state, but a resonance in the 27r continuum. In order to determine whether the continuum 
contribution in the instanton model is predominantly 27r or 2 quark would require the determination 
of the corresponding three point functions (which has not been done yet). Instead, we will compare the 
full correlation function with the non-interacting one, where we use the average (constituent quark) 
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propagator determined in the same ensemble (Fig. As explained above, this comparison provides 
a measure of the strength of interaction. We observe that there is an attractive interaction generated 
in the interacting liquid. The interaction is due to correlated instanton-anti-instanton pairs. This is 
consistent with the fact that the interaction is considerably smaller in the random ensemble. In the 
random model, the strength of the interaction grows as the ensemble becomes more dense. However, 
the interaction in the full ensemble is significantly larger than in the random model at the same 
diluteness. Therefore, most of the interaction comes from dynamically generated pairs. 

The situation is drastically different in the r/ channel. Among the ~ 40 correlation functions 
calculated in the random ensemble, only the Vj and the isovector-scalar 8 were found to be completely 
unacceptable: The correlation function decreases very rapidly and becomes negative at x ~ 0.4 fm. 
This behavior is incompatible even with a normal spectral representation. The interaction in the 
random ensemble is too repulsive, and the model "over-explains" the \J(X)a anomaly. 

The results in the unquenched ensembles (closed and open points) significantly improve the sit- 
uation. This is related to dynamical correlations between instantons and anti-instantons (topological 
charge screening). The single instanton contribution is repulsive, but the contribution from pairs is 
attractive. Only if correlations among instantons and anti-instantons are sufficiently strong, the cor- 
relators are prevented from becoming negative. Quantitatively, the (5 and 77ns masses in the streamline 
ensemble are still too heavy as compared to their experimental values. In the ratio ansatz, on the 
other hand, the correlation functions even show an enhancement at distances on the order of 1 fm, 
and the fitted masses are too light. This shows that the r/ channel is very sensitive to the strength of 
correlations among instantons. 

In summary, pion properties are mostly sensitive to global properties of the instanton ensemble, 
in particular its diluteness. Good phenomenology demands p^n ~ 0.03, as originally suggested in.d 
The properties of the p meson are essentially independent of the diluteness, but show sensitivity to 
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Fig. 6. Rho meson correlation functions. The various curves and data sets are labeled as in Fig. ^ The dashed 
squares show the non-interacting part of the rho meson correlator in the interacting ensemble. 
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// correlations. These correlations become crucial in the r( channel. 
5.^. Baryonic correlation functions 

Existence of a strongly attractive interaction in the pseudoscalar quark-antiquark (pion) channel 
also implies an attractive interaction in the scalar quark-quark (diquark) channel. This interaction is 
phenomenologically very desirable, because it immediately explains why the nucleon and lambda are 
light, while the delta and sigma are heavy. 

The proton currents (with nOpderivatives and the minimum number of quark fields) with positive 
parity and spin 1/2 are given byBj It is useful to rewrite these currents in terms of scalar and 
pseudoscalar diquarks 

771,2 = (2,4){e„b,(u'^Crf^)75?/=Tea6cKC75d'')u^}. (37) 

Nucleon correlation functions are defined by 11^^ (x) = (r/Q(0)fy^(a;)), where a, (3 are the Dirac indices 
of the nucleon currents. In total, there are six different nucleon correlators: the diagonal 771771, 772772 
and off-diagonal 771772 correlators, each contracted with either the identity or 7 • a;. Let us focus on the 
first two of these correlation functions (for more detail, seeQ and references therein). 





Fig. 7. Nucleon and delta correlation functions 11^ and 11^. Curves labeled as in Figs. 

The correlation function 11^ in the interacting ensemble is shown in Fig. 0. There is a significant 
enhancement over the perturbative contribution which is nicely described in terms of the nucleon 
contribution. Numerically, we find[J] tun — 1.019 GeV. In the random ensemble, we have measured the 
nucleon mass at smaller quark masses and found = 0.96 ± 0.03 GeV. The nucleon mass is fairly 
insensitive to the instanton ensemble. However, the strength of the correlation function depends on the 
instanton ensemble. This is reflected by the value of the nucleon coupling constant, which is smaller 

tNote that this value corresponds to a relatively large current quark mass m = 30 MeV. 



18 



in the IILM. InE^ we studied all six nucleoli correlation functions. We showed that all correlation 
functions can be described with the same nucleon mass and coupling constants. 

The fitted value of the threshold is Eq ~ 1.8 GeV, indicating that there is little strength in the 
"three quark continuum" (dual to higher resonances in the nucleon channel) . The fact that the nucleon 
in IILM is actually bound can also be demonstrated by comparing the full nucleon correlation function 
with that of three non-interacting quarks (the cube of the average propagator. The full correlator is 
significantly larger than the non-interacting one. 

Significant part of this interaction was traced down to strongly attractive scalar diquarks. The 
nucleon (at least in IILM) is a strongly bound diquark, plus a loosely bound the third quark. 

The properties of this diquark picture of the nucleon continue to be disputed by phenomenol- 
ogists. We would return to diquarks in the next section, where they would become Cooper pairs of 
Color Superconductors. 

In the case of the A resonance, there exists only one independent current, given (for the A"''+) 

by 

= eabc{u''C-fy)u^. (38) 

However, the spin structure of the correlator n^j^.^^(a;) — (?7^a (0)77^(2;)) is much richer. In general, 
there are ten independent tensor structures, but the Rarita-Schwinger constraint 7^77^ — reduces 
this number to four. 

The mass of the delta resonance is too large in the random model, but closer to experiment 
in the unquenched ensemble. Note that similar to the nucleon, part of this discrepancy is due to the 
value of the current mass. Nevertheless, the delta-nucleon mass splitting in the unquenched ensemble 
is TOA — AT = 409 MeV, larger but comparable to the experimental value 297 MeV. It mostly comes 
from the absence of attractive scalar diquarks in this channel. 
6. The Phases of QCD 

Turning to finite temperature/density QCD now, let me start with emphasizing its main goals. 
Those are not to use once again a semi-classical or perturbative calculations, similar to what have 
been done before in vacuum. What we are looking for here are new phases of QCD (and related 
theories), namely new self-consistent solutions which differs qualitatively from what we have in the 
QCD vacuum. 

One such phase occurs at high enough temperature T > it is known as Quark Glumi Plasma 
(QGP). It is a phase understandable in terms of basic quark and gluon-like excitations,Efl without 

confinement and with unbroken chiral symmetry in the massless limitf]. One of the main goals of 
heavy ion program, especially at new Brookhaven dedicated facility RHIC, is to study transitions to 
this phase. 

Another one, which gets much attention recently, is the direction of finite density. Very robust 
Color Superconductivity was found to be the case here. Let me also mention one more frontier not 
to be discussed below, which has not yet attracted sufficient attention: namely transition (or many 
transitions?) as the number of light flavors Nf grows. The minimal scenario includes transition from 
the usual hadronic phase to one more unusual QCD phase, the conformal one, in which there are no 
particle-like excitations and correlators are power-like in the infrared. Even the position of the critical 
point is unknown. 

That the main driving force of these studies is the intellectual challenge it provides. In order to 
prove to you that the problems related with mechanism of conflnement and chiral symmetry breaking 
are not simple, it is enough to remind that a lot of people who moved these days into quantum 
gravity /string theory did so partly because of frustration with the non-perturbative QCD. However I 
still think that, with the help of all experimental data, real and numerical, we have at hand for QCD 
we still have much better chances here. 

t It docs not mean though, that it is a simple issue to understand even the high-T limit of QCD, related to non- 
perturbative 3d dynamics. 
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6.1. The Phase diagram 

The QCD phase diagram version circa 19990 is shown in Fig ||(a), at the baryonic chemical 
potential /x (normalized per quark, not per baryon) - the temperature T plane. Some part of it is old: 
it is hadronic phase at small values of both parameters, and QGP phase at large T,/i. 

The phase transition line separating them most probably does not really start at T = Tc, A* = 
but at a "endpoint" E, a remnant of the so called QCD tricritical point which QCD has in the chiral 
(all quarkS|-aj"e massless) limit. Although we do not know where it isQ, we hope we know how to 
find it, see.EJ The proposed ideas rotate around the fact that the order parameter, the VEV of the 
sigma meson, is at this point truly massless, and creates a kind of a "critical opalesence" . Similar 
phenomena were predicted and then indeed observed at the endpoint of another line (called M from 
multi- fragmentation), separating liquid nuclear matter from nuclear gas phase. 

The large-density (and low-T) region looks rather different from what was shown at conferences 
just few months ago: there appear two new Color Super-conducting phases there. Unfortunately heavy 
ion collisions do not cross this part of the phase diagrams, and so it belongs to a neutron star physics. 




Fig. 8. (a) Schematic phase diagram of QCD, in temperature T- baryon chemical potential fi plane. E and M 
show critical endpoints of first order transitions: M (from multi-fragmentation) is that for liquid-gas transition 
in nuclear matter. The color superconducting phases, CSC2 and CSC3 are explained in the text, (b) The 
part of the phase diagram which has been experimentally accessible by heavy ion collisions. Closed (open) 
points correspond to chemical (thermal) freeze-out. The zigzag is the adiabatic paths correspond to entropy 
per baryon ratio indicated. Lines of thermal freeze-out for central sulfur-sulfur and PbPb collisions, as well as 
the initial line for PbPb are also indicated. 

The 2-flavor-like color supercpiuluctor CSC2 phase was known before,0 but realization that 
it should be induced by instantonsc3cil has increased the gaps (and Tc) from a few MeV scale to 
~100 MeV (50 MeV). In the hindsight is should hardly be surprising, since the same interaction 
in qq channel is responsible for chiral symmetry breaking, producing the gap (the constituent quark 
mass) as large as 350-400 MeV. Furthermore, in the two-color QCD, there is the so called Pauli-Gursey 
symmetry which relates these two condensates. So at high density the chiral condensate < qq > simply 
rotates into superconductor one < qq >, while the gap remains the same. 

The symmetries of the CSC2 phase are similar to the electroweak part of the Standard Model, 
with the condensed scalar isoscalar ud diquark operating as Higgs. The colored condensate breaks 
the color group, making 5 out of 8 gluons massive. The 3-flavor-like phase, CSC3, is brand new: it 

§Its position is very sensitive to precise value of the strange quark mass nis 
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was proposed iiJl3 based on one gluon exchange interaction, but in fact it is favored by instantons as 
wello Its unusual features include color-flavor locking and coexistence of both types of condensates, 
< qq > and < qq >. It combines features of the Higgs phase (8 massive gluons) and of the usual 
hadronic phase (8 massless "pions"). 

Above I mention approach to high density starting from the vacuum. One can also work out m 
the opposite direction, starting from very large densities and going down. As it was noticed in ref.,E3 
electric part of one-gluon exchange is screened, and therefore the Cooper pairs appears due to magnetic 
forces. It is interesting by itself, as a rare example: one has to take care of time delay effects of the 
interaction. The result is the indefinitely growing gaps at large /i > \QGeV , as A ~ ixexp{— -^-^) ■ 

6. 2. Finite T transition and Large Number of Flavors 

There is no place here to discuss this subject in details: there are rather extensive lattice data 
now, and so we actually know quite a lot about these transition. 

Let me only emphasize what it is look like from the perspective of the instanton-based theory. If 
the near-random set of instantons leads to chiral symmetry breaking and quasi-zero modes at low T, 
we should be able to explain in the same terms how the high-T phase look like. The simplest solution 
would be just disappearance of instantons at T > T^, and at some early time people thought this 
is what actually happens. However, it should not be like this because the Debye screening which is 
killing them only appears atT — Tc- Lattice data works have also found no depletion of the instanton 
density up to T = Tc. 

On the other hand, the absence of the condensate and quasi-zero modes can only mean that 
the "liquid" is now broken into finite pieces. The simplest of them are pairs, or the instanton-anti- 
instanton molecules. This is precisely what instanton simulations have found.Q Whether it is or is not 
so on the lattice is not yet clear. Some nice molecules were located, but the evidences for the molecular 
mechanism are still far from being completely convincing. No alternative have been so far proposed, 
however. 

The results of simulations with Nf — 2, 3, fiavors with equal masses can be summarized 
as follows. For Nf = 2 there is second order phase transition which turns into a line of first order 
transitions in the m — T plane for Nf > 2. If the system is in the chirally restored phase (T > Tc) 
at m = 0, we find a discontinuity in the chiral order parameter if the mass is increased beyond some 
critical value. Qualitatively, the reason for this behavior is clear. While increasing the temperature 
increases the role of correlations caused by fermion determinant, increasing the quark mass has the 
opposite effect. We also observe that increasing the number of fiavors lowers the transition temperature. 
Again, increasing the number of fiavors means that the determinant is raised to a higher power, so 
fermion induced correlations become stronger. For Nf = 5 we find that the transition temperature 
drops to zero and the instanton liquid has a chirally symmetric ground state, provided the dynamical 
quark mass is less than some critical value. Studying the instanton ensemble in more detail shows that 
in this case, all instantons are bound into molecules. 

Unfortunately, little is known about QCD with different numbers of flavors from lattice simu- 
lations. There are data by the Columbia group for Nf = 4. The most important result is that chiral 
symmetry breaking effects were found to be drastically smaller as compared to Nf — 0,2. In particu- 
lar, the mass splittings between chiral partners such as tt — cr, p ~ ai, N{^^) — ); extrapolated 
to m = were found to be 4-5 times smaller. This agrees well with what was found in the interacting 
instanton model: more work in this direction is certainly needed. 

6.3. Color superconductivity 

My interest was initiated by findinjsl that in the instanton liquid model even without any 
quark matter, the ud scalar diquarks are very deeply bound, by amount comparable to the constituent 

case N f = 4 is omitted because in this case it is very hard to determine whether the phase transition happens at 

T > 0. 
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quark mass. So, phenomenological manifestationsEj of such diquarks have in fact deep dynamical 
roots: they follow from the same basic dynamics as the "superconductivity" of the QCD vacuum, the 
chiral (x-)symmetry breaking. These spin-isospin-zero diquarks are related to pions, and should be 

quite robust element of nucleon (octet baryons) structure|||] . 

Another argument for deeply bound diquarks comes from bi-color {Nc — 2) theory: in it the 
scalar diquark is degenerate with pions. By continuity from Nc = 2 to 3, a trace of it should exist in 
real QCDQ. 

Explicit calculations with instanton-induced forces for Nf — 2,Nc = 3 QCD have been made 

in two simultaneous^ papers. Indeed, a very robust Cooper pairs and gaps A ~ 100 MeV were 
found. From then on, the field is booming. 

Instantons create the following amusing triality: there are three attractive channels which com- 
pete: (i) the instanton-induced attraction in qq channel leading to ^-symmetry breaking, (ii) the 
instanton-induced attraction in qq which leads to color superconductivity, (iii) the light- quark-induced 
attraction of II, which leads to pairing of instantons into "molecules" and a Quark-Gluon Plasma 
(QGP) phase without any condensates. 

How the calculations are actually made ?. Analytically, mostly in the mean field approximation, 
similar to the original BCS theory in Gorkov formulation. Total thermodynamical potential consists 
of "kinetic energy" of the quark Fermi gas, including mass operators of two types (shown in figure 
below). The "potential energy" in such approximation is the interaction Lagrangian convoluted with 
all possible condensates. For example, instanton-induced one with Nf = 3 leads to two types of 
diagrams shown in Fig. 4, with (a) < qq >^ and (b) < qq >^< qq >. Then one minimizes the potential 
over all condensates and get gap equations: algebra may be involved because masses/condensates are 
color-flavor matrices. 




Fig. 9. (a) Two mass operators; (b) and example of two contribution to the potential energy for 3 flavor theory, 
including cube of the < qq > and | < qq > \'^ < qq >. 



Bi-color QCD: a very special theory One reason it is special is well known to lattice 
community: its fermionic determinant is real even for non-zero fj,, which makes simulations possible. 

However the major interest to this theory is related the so called Pauli-Gursey symmetry. We 
have argued above that pions and diquarks appear at the same one-instanton level, and are so to say 
brothers. In bi-color QCD they becomes identical twins: due to additional symmetry mentioned the 
diquarks are degenerate with mesons. 

In particular, the x-symmetry breaking is done like this SU{2Nf) Sp{2Nf), and for Nf — 2 
the coset K — SU{4)/Sp{4) = SO{6)/SO{5) ~ S^. Those 5 massless modes are pions plus scalar 
diquark S and its anti-particle S. 

Vector diquarks are degenerates with vector mesons, etc. Therefore, the scalar-vector splitting 
is in this case about twice the constituent quark mass, or about 800 MeV. It should be compared to 
binding in the "real" Nc = 3 QCD of only 200-300 MeV, and to zero binding in the large-iVc limit. 

II As opposed to A (decuplet) baryons. 

** Instanton-induced interaction strength in diquark channel is l/{Nc — 1) of that for q-f^q one. It is the same at Nc = 2, 
zero for large N^, and is exactly in between for Nc = 3. 
Submitted to hep-ph on the same day. 
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The corresponding_aigma model describing this x-symmetry breaking was worked out in@: for 
further development see.E3 As argued in^J in this theory the critical value of transition to Color 
Superconductivity is simply /i ~ 711^^/2, or zero in the chiral limit. The diquark candensate is just 
rotated < qq > cpe, and the gap is the constituent quark mass. Recent lattice workalj and instanton 
liquid simulationEil display it in great details, building confidence for other cases. 

Two flavor Three coloi;i-.QCp: the CSC2 phase The first studies of instanton-induced 
CSC were made for this theory :.E2lEilEa The gap around fOO MeV was found, and the main point is 
a competition between the usual x-symmetry breaking (pairing in the qq channel) with CSC (pairing 
in the qq channel). In the corresponding high density phase we called CSC2 the chiral symmetry is 
simply restored. 

In all these works one more possible phase (intermediate between vacuum and CSC2), Fermi gas. 
of constituent quarks, with both M, A ^ - was unstable. However in last more refined calculationea 
it obtains a small window. Its features are amusingly close to those of nuclear matter: but it. isn't, of 
course: to get nucleons one should go outside the mean field. First attempted to do so ircJ was for 
another cluster - the // molecules. At T=0 it is however only 10% correction to previous results, but 
is dominant as T grows. _ 

Three color Three flavor QCD produces a new phenomenon called color-flavor lockings^ It 
means that the diquark condensate has the following structure (qfCqj) = AiSiaSbj + A25ibSja, where 
ij are color and ab fiavor indices. It is very symmetric, reducing SU{3)cSU{3) f — * SU (3) diagonal- We 
call it the CSC3 phase: it still breaks chkal symmetry. 

It was verified to be the lowest iriEj for the one gluon exchange interaction (adequate for high 
density): in this case only jqq^ condensate exists. About .the same phase appears for the instanton- 
induced Lagrangian (adequate for intermediate density)Ej: but now the < qq > condensates are 
non-zero as well. The color-flavor locking is probably always the case for that theory. 

Gaps Si and masses ai (proportional to < gg >), following from instanton-based calculation,Cj 
are shown as a function of fj, in the following figure 
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Physics issues under discussion for CSC3 include hadron-quark continuity. As pointed out in,E3 
the CSC3 phase not only has the same symmetries as hadronic matter (e.g. broken x-synimetry) , 
but also very similar excitations. 8 gluons become 8 massive vector mesons, 3*3 quarks become 8-1-1 
"baryons". The 8 massless pions remain masslessf^ Furthermore, photon and gluons are combined 

MVery exotic 3d objects, "super-qualitons" ,t3 the skyrmions made of pions are among the excitations. 
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into a massless ^inside- Calculation of masses and coupling constants of all of them is now in progress. 
Can these phases be distiquished, and should there be any phase transition in the Nf — 3 theory, 
separating it from nuclear matter? There is no need for it, at least from symmetry point nf view 

The realistic QCD(Two light plus strange flavor) was studied in several papers.aEj Just 
kinematically it is easy to see that us,ds Cooper pairs with zero momentum is difhcult to make: 
for fiu,d — IJ's the momenta ^ pf . Instantons generate also a dynamical operator ms{ud){ud). 
Resulting behavior is as shown in our first figure. l-l 

Asymptotically large densities At /i >> IGeV the instantons are Debye-screenedjil as well 
as the electric (Coulomb) gluons. So magnetic gluons overtake electric ones.cZl Magnetically bound 
Cooper pair is interesting by itself, as a rare example: one has to take care of time delay effects with 
Eliashberg eqn, etc. Angular integral leads to second log in the gap equation, leading to unusual 
answer: A ~ /i exp{~3Tr'^ / \/2g) which implies that the gap grows indefinitely with /i Q and pQCD 
becomes finally justified. However, it is the case for huge densities, with fj, > lOGeV or so. 

Final remarks. All variants of CSC have an "internal photon" (a combination of the photon 
and gluon) for which the condensate is not charged and which therefore is not expelled from it. So, is 
it a superconductor, after all? 

I think the answers is still "yes". For example, if one puts a piece of CSC into a magnet, it 
should still levitate: although ^inside is massless, the magnet uses ^outside held and a part of it is 
expelled^ 

Let me finish this section with few homework questions. What is the role of confinement in 
all these transitions? What is nuclear matter for different quark masses, anyway? Do we have other 
phases in between, like diquark-quark phase or (analog of) K condensation, or different crystal-like 
phases? Is there indeed a (remnant of) the tricritical point which we can find experimentally? And, 
How can we do finite density calculations on the lattice? 
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